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Abstract. These notes are devoted to the statistical mechanics of directed polymers 
interacting with one-dimensional spatial defects. We are interested in particular in the 
situation where frozen disorder is present. These polymer models undergo a localiza- 
tion /derealization transition. There is a large (bio)-physics literature on the subject 
since these systems describe, for instance, the statistics of thermally created loops in 
DNA double strands and the interaction between (1 + l)-dimensional interfaces and 
disordered walls. In these cases the transition corresponds, respectively, to the DNA de- 
naturation transition and to the wetting transition. More abstractly, one may see these 
models as random and inhomogeneous perturbations of renewal processes. 

The last few years have witnessed a great progress in the mathematical understanding 
of the equilibrium properties of these systems. In particular, many rigorous results about 
the location of the critical point, about critical exponents and path properties of the 
polymer in the two thermodynamic phases (localized and delocalized) are now available. 

Here, we will focus on some aspects of this topic - in particular, on the non-perturbative 
effects of disorder. The mathematical tools employed range from renewal theory to large 
deviations and, interestingly, show tight connections with techniques developed recently 
in the mathematical study of mean field spin glasses. 
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1. Introduction and motivations 



Consider a Markov chain {SnjneN on some state space Q, say, $7 = Z . We can unfold S 
along the discrete time axis, i.e., we can consider the sequence {(n, SVOImeN and interpret 
it as the configuration of a directed polymer in the space N x J7. In the examples which 
motivate our analysis, the discrete time is actually better interpreted as one of the space 
coordinates. The "directed" character of this polymer just refers to the fact that the 
first coordinate, n, is always increasing. In particular, the polymer can have no self- 
intersections. Some assumptions on the law of the Markov chain will be made in Section 
[21 where the model is defined precisely. Now let be a specific point in f2, and assume that 
the polymer receives a reward e (or a penalty, if e < 0) whenever S n = 0, i.e., whenever 
it touches the defect line N x {0}. In other words, the probability of a configuration of 
{S\, S2, • • • , Sn} is modified by an exponential, Boltzmann-type factor 



It is clear that if e > contacts with the defect line are enhanced with respect to the 
e = (or free) case, and that the opposite is true for e < 0. One can intuitively expect 
that the in the thermodynamic limit N — > 00 a phase transition occurs: for e > e c the 
polymer stays close to the defect line essentially for every n, while for e < e c it is repelled 
by it and touches it only at a few places. This is indeed roughly speaking what happens, 
and the transition is given the name of localization/ 'derealization transition. We warn 
the reader that it is not true in general that the critical value is e c = 0: if the Markov 
chain is transient, then e c > 0, i.e., a strictly positive reward is needed to pin the polymer 
to the defect line (cf . Section I2.6H . 

A more interesting situation is that where the constant repulsion/attraction e is replaced 
by a local, site-dependent repulsion/attraction e n . One can for instance consider the 
situation where e n varies periodically in n, but we will rather concentrate on the case 
where e n are independent and identically distributed (IID) random variables. We will 
see that, again, the transition exists when, say, the average e of e n is varied. However, 
in this case the mechanism is much more subtle. This is reflected for instance in the 
counter-intuitive fact that e c may be negative: a globally repulsive defect line can attract 
the polymer! Presence of disorder opens the way to a large number of exciting questions, 
among which we will roughly speaking select the following one: how are the critical point 
and the critical exponents influenced by disorder? 

There are several reasons to study disordered pinning models: 

• there is a vast physics and bio-physics literature on the subject, with intriguing 
(but often contradictory) theoretical predictions and numerical/experimental ob- 
servations. See also Section \TQ\ 

• they are interesting generalizations of classical renewal sequences. From this point 
of view they raise new questions and challenges, like the problem of the speed of 
convergence to equilibrium for the renewal probability in absence of translation 
invariance (cf. in particular Section [6]) ; 

• finally (and this is my main motivation) they are genuinely quenched-disordered 
systems where randomness has deep, non-perturbative effects. With respect to 
other systems like disordered ferromagnets or spin glasses, moreover, disordered 
pinning models have the advantage that their homogeneous counterparts are under 
full mathematical control. These models, therefore, turn out to be an ideal testing 
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ground for theoretical physics arguments like the Harris criterion and renormal- 
ization group analysis. 

It is also quite encouraging, from the point of view of mathematical physics, that rigorous 
methods have been able not only to confirm predictions made by theoretical physicists, 
but in some cases also to resolve controversies (it is the case for instance of the results in 
Section \5.6\ which disprove some claims appeared previously in the physical literature). 

1.1. A side remark on literature and on the scope of these notes. A excellent 
recent introductory work on pinning models with quenched disorder (among other topics) 
is the book [22\ by Giambattista Giacomin. In order to avoid the risk of producing a 
resume of it, we have focussed on aspects which are not (or are only tangentially) touched 
in [22]. On the other hand, we will say very little about "polymer path properties", 
to which Chapters 7 and 8 of [22] are devoted. A certain degree of overlap is however 
inevitable, especially in the introductory sections [2] and [U results taken from [22] will be 
often stated without proofs (unless they are essential in the logic of these notes). 

We would also like to mention that some of the results of these notes apply also to a 
model much related to disordered pinning, namely random heteropolymers (or copolymers) 
at selective interfaces. It is the case, for instance, of the results of Sections 15.61 and l6l We 
have chosen to deal only with the pinning model for compactness of presentation, but we 
invite readers interested in the heteropolymer problem to look, for instance, at [11] . [37] . 
[22] and references therein. 



2. The model and its free energy 

2.1. The basic renewal process ("the free polymer"). Our starting point will be 
a renewal r on the integers, r := {Tj}j = o,i,2,...> where tq = and {ti — Tj_i}j>i are IID 
positive and integer-valued random variables. The law of the renewal will be denoted by 
P, and the corresponding expectation by E. In terms of the "directed polymer picture" 
of the introduction, P is the law of the set r of the points where the polymer touches the 
defect line, in absence of interaction: r = {n : S n = 0} (cf. also Section [2U|) . We assume 
that (Tj — Tj_i) or, equivalently, t\ is P-almost surely finite: if 

Kin) := P(n = n), (2.1) 

this amounts to requiring ]CneN-^( n ) = This, of course, implies that the renewal is 
recurrent: P-almost surely, r contains infinitely many points. A second assumption is 
that K(.) has a power-like tail. More precisely, we require that 

K{n) = for every n G N, (2.2) 

for some a > and a slowly varying function L(.). We recall that a function (0,oo) 9 
x — ► L(x) £ (0, oo) is said to be slowly varying at infinity if [8] 

Km ^ = 1 (2.3) 

- L(x) 



- oc 



for every r > 0. In particular, a slowly varying function diverges or vanishes at infinity 
slower than any power. The interested reader may look at [8j for properties and many 
interesting applications of slow variation. Of course, every positive function L(.) having a 
non-zero limit at infinity is slowly varying. Less trivial examples are L(x) = (log(l + x)) 1 
for 7 G R. 
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Observe that the normalization condition ^2 ne jqK(n) = 1 implies that, if a = 0, L(.) 
must tend to zero at infinity (cf. also Section 12.61 below for an example) . 

It is important to remark that typical configurations of r are very different according to 
whether a is larger or smaller than 1. Indeed the average distance between two successive 
points, 

B(T i -T i _ 1 ) = J2nK(n), (2.4) 

is finite for a > 1 and infinite for a < 1. In standard terminology, r is positively recurrent 
(i.e., t occupies a finite fraction of N) for a > 1 and null-recurrent for a < 1 (the density 
of r in N is zero). This is a simple consequence of the classical renewal theorem [61 Chap. 
I, Th. 2.2], which states that 

The distinction a ^ 1 plays an important role, especially in the behavior of the homoge- 
neous pinning model (cf. Section [3]). Later on we will see the emergence of an even more 
important threshold value: a c = 1/2. 

Remark 2.1. For a = 1, the question whether the renewal is positively or null recurrent is 
determined by the behavior at infinity of L(.): from (|2.5j) we see that t is finitely recurrent 
iff X^n-^( n )/ n < 00 • F° r i ns t ance > one has null recurrence if L{.) has a positive limit at 
infinity. 

2.2. The model in presence of interaction. Now we want to introduce an interaction 
which favors the occurrence of a renewal at some points and inhibits it at others. To this 
purpose, let u> (referred to as quenched randomness or random charges) be a sequence 
{uJ n } n <=N of IID random variables with law P. The basic assumption on u n , apart from 
the fact of being IID, is that Ewi = and Ew? = 1. These are rather conventions 
than assumptions, since by varying the parameters (5 and h in Eq. (12. 6p below one can 
effectively tune average and variance of the charges. To be specific, in these notes we will 

consider only two (important) examples: the Gaussian case u\ = J\f(0, 1) and the bounded 
case, \ui\ < C < oo. Many results are expected (or proven) to hold in wider generality 
and a few remarks in this direction are scattered throughout the notes. 

We are now ready to define the free energy of our model: given h 6 R, (3 > and A'Gff 

let 

*%(P,h) ■= ^log Z N ,Mh) := llogE (e^^ +h ^5 N ) , (2.6) 

where for notational simplicity we put 5 n := l{ n g T }, 1a being the indicator function of 
a set A. The quenched average of the free energy, or quenched free energy for short, is 
defined as 

F N (J3,h):=EF%(p,h). (2.7) 

Note that the factor 8n in (|2.6p corresponds to imposing the boundary condition N £ r 
(the boundary condition G r at the left border is implicit in the law P). One could 
equivalently work with free boundary conditions at N (i.e., replace 5n by 1). The infinite- 
volume free energy would not change, but some technical steps in the proofs of some results 
would be (slightly) more involved. 
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We need also a notation for the Boltzmann-Gibbs average: given a realization u of the 
randomness and a system size N, for a P-measurable function /(.) set 

E (7(r) eESLiOa^+fc)*™^ 



2.3. Existence and non-negativity of the free energy. As usual in statistical me- 
chanics, one is (mostly) interested in the thermodynamic limit (i.e., the limit ./V — * oo). 
A classical question concerns the existence of the thermodynamic limit of the free energy, 
and its dependence on the realization of the randomness u. In the context of the models 
we are considering, the answer is well established: 

Theorem 2.2. [22l Th. 4.1] 7/E|wi| < oo, the limit 

F(p,h):= lim ±-log Z N p(J3,h) (2.9) 

N—+co iv 

exists for every (3 > 0, h G M and it is P( dto)-almost surely independent of u. 

Of course, the limit does depend in general on the law P of the disorder. 

Note that the only assumption on disorder, apart from the IID character of the charges, 
is finiteness of the first moment, so that existence and self- averaging of the infinite- volume 
free energy holds in much wider generality than in the cases of Gaussian or bounded 
disorder we are considering here. 

Some properties of the free energy come essentially for free: in particular, F({3, h) 
is convex in ((3,h), non-decreasing in h, continuous everywhere and differentiable al- 
most everywhere as a consequence of convexity. Another easy fact is that the sequence 
{N F N (J3, h)} N( z N is super-additive: for every N,M e N, one has (N + M)F N+M (P, h) > 
NFn(J3, h) + MFm(P, h). This is easily proven: 

(N + M)F N+M (p,h) = ElogE (e^ U ^ +h ^S N+M ) (2.10) 

> E log E (e^^ +h ^S N ^^ + ^ +h ^S N+M ) 
= NF N ((3,h) + MF M ((3,h), 

where in the last step we used invariance of P with respect to left shifts and the renewal 
property of P. It is a standard fact that super-additivity implies 

F(p, h) > F N {p, h) for every G N. (2.11) 

2.4. Contact fraction and critical point. As we already mentioned, the interest in 
this class of models is mainly due to the fact that they show a so-called localization- 
delocalization transition. This is best understood in view of the elementary bound F(/3, h) > 
0. This positivity property is immediate to prove: 

F N ((3,h) > ^ElogE (e^i(^ + ^l {T1=iV} ) = A + Ll og K(N) (2.12) 

and the claimed non-negativity in the limit follows from (|2.2p . Recalling that F(/3,h) is 
non-decreasing in h, for a given (5 the localization/delocalization critical point is defined 
to be 

hc{P) ■= sup{h : F{p, h) = 0} (2.13) 
and the function (3 — ► h c {(5) is referred to as the critical line. The region of parameters 

C:={(J3,h):p>0,h>h c (J3)} 
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and 

V:={{P,h):0>O,h<h e (0)} 

are referred to as localized and delocalized phases, respectively. Since level sets of a 
convex function are convex, £ is a convex set and the function h c {.) : [0, oo) 3 (3 — » h c {(5) 
is concave. The reason for the names "localized" and "delocalized" can be understood 
looking at the so-called contact fraction £n, defined through 

fa:= |m{i...,* } | (2 , 4) 

and taking values between and 1 (as usual, \A\ denotes the cardinality of a set A). It is 
immediate to check that 

d h F%([3,h)=-E^(£ N ) (2.15) 

and, by standard arguments based on convexity, this equality survives in the thermody- 
namic limit whenever the free energy is differentiable: 

lim E JI(£n) =' d h F((3, h) for every h such that dtF((3, h) = d7F((3, h). (2.16) 

N—>oo ' 

We have already mentioned that differentiability holds for Lebesgue-almost every value 
of h. However, much more than this is true: as it was proven in [25J, differentiability 
(actually, infinite differentiability) in h holds whenever h > h c (/3). We can therefore 
conclude the following: for h < h c (j3) (or for h < h c {(3) if F(f3,h) is differentiable at 
h c (P)) the thermal average of the contact fraction tends for to zero for N — > oo (almost 
surely in the disorder), while for h > h c (f3) it tends to dhF((3, h) > 0. The average contact 
fraction plays the role of an order parameter, like the spontaneous magnetization in the 
Ising model, which is zero above the critical temperature and positive below it. 

Actually, much more refined statements about the behavior of the contact fraction in 
the two phases are available. In particular: 

• for statements concerning the localized phase we refer to [25] . There, it is proven 
that, roughly speaking, not only typical configurations r have a number 

N£ N ~Nd h F(0,h) 

of points, but also that these points are rather uniformly distributed in {1, ... , N}: 
long gaps between them are exponentially suppressed, and the largest gap is of 
order logiV (cf. Theorem 16.31 below) : 

• for h < h c ({3) we refer to [23] and [22j Ch. 8], where it is proven that £n is typically 
at most of order (log N)/N. 

In this sense, if one goes back to the pictorial image of r as the set of points of polymer- 
defect contact, one sees that the definition of (de)-localization in terms of free energy, as 
given above, does indeed correspond to the intuitive idea in terms of path properties: in 
C the polymer stays at distance O(l) from the defect, while in D it wanders away from it 
and touches it only a small (at most log iV) number of times. 

The reader should remark that we have made no conclusive statement about the be- 
havior of the contact fraction at h c {(5), since we have not attacked yet the very important 
question of the regularity of the free energy at the critical point. This will be the subject 
of Sections HI and El 
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2.5. Quenched versus annealed free energy. Inequality (|2.12|) is a good example of 
how selecting a particular subset of configurations (in that case, those for which t\ = N) 
provides useful free energy lower bounds. For more refined results in this direction we 
refer to [5] and [22j Sec. 5.2]. There, this technique is employed to prove that h c (j3) is 
strictly decreasing as a function of (3 which implies in particular that, since h c {.) is concave, 
h c {(3) tends to — oo for (3 — ► oo. This corresponds to the apriori non-intuitive fact that, as 
mentioned in the introduction, even if the charges are on average repulsive the defect line 
can pin the polymer. This is purely an effect of spatial inhomogeneities due to disorder: 
for (3 large, it is convenient for the polymer to touch the defect line in correspondence of 
attractive charges, where it gets a reward f3iv n + h >> 1, while the entropic cost of avoiding 
the repulsive charges is independent of [3. Free energy lower bounds were obtained also in 
the study of a different model, the heteropolymer at a selective interface, in [10J. 

Free energy upper bounds are on the other hand more subtle to get. An immediate one 
can be however obtained by a simple application of Jensen's inequality: 

F N (/3,h) <±logEZ NjU1 (f3,h) = ^logE (eZZ =1 (h+io g M m Sn §N ^ (217) 

= F N (0,h + ]ogM(P)) =:Fft(p,h), 

where M(/3) := Ee^ 1 . In particular, log M(f3) = f3 2 /2 in the case of Gaussian disorder. 
F a (f3, h) := F(0, (3 + log M{(3)) is referred to as annealed free energy, and we see that it is 
just the free energy of the homogeneous system (with the same choice of K(.)) computed 
for a shifted value of h. The physical interpretation of the annealed free energy is clear: 
since configurations of to and r are averaged on the same footing, it corresponds to a 
system where impurities can thermalize on the same time-scales as the "polymer degrees 
of freedom" (i.e., r). This is not the physical situation one wishes to study (quenched 
disorder corresponds rather to impurities which are frozen, or which can evolve only on 
time-scales which are so long that they can be considered as infinite from the experimental 
point of view). All the same, the information provided by (|2.17p is not at all empty. Define 
first of all the annealed critical point as 

h*{fi) := snp{h : F a ((3, h) = 0} = h c (0) - log M(/3). (2.18) 

Thanks to (12.171) and (12. 13ft . one has immediately 

h c ((3) > h c (0)- log M(f3), (2.19) 

a bound which, as will be discussed in Section [5.31 is optimal for a < 1/2 and (3 small. 

2.6. Back to examples and motivations. Typical examples of renewal sequences sat- 
isfying (|2.ip . (|2.2p are the following. Let {S n } n >o be the simple random walk (SRW) on 
Z, with law P SRW and Sq := 0, i.e., {S n — S n -\} n ^ are HD symmetric random vari- 
ables with values in { — 1, +1}- Then, it is known that |19j t := {n £ N : S^n = 0} is a 
null-recurrent renewal sequence such that the law of t\ satisfies (|2.2p with a = 1/2 and 
L(.) asymptotically constant. The reason why one looks only at even values of n in the 
definition of r in this case is due just to the periodicity of the SRW. If instead one takes 
the SRW on Z, 2 , then r (defined exactly as above) is always a null-recurrent renewal but 
in this case a = and L(n) ~ c/(logn) 2 [33]. Note that in this case, the presence of the 
slowly varying function L(.) is essential in making K (.) summable. 

What happens in the case of the SRW on Z d when d > 3? This example does not 
fall directly into the class we are considering since this process is transient, and therefore 
the set r of its returns to zero is a transient renewal sequence. However this is not too 
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bad. Indeed, suppose more generally that one is given K(.) which satisfies (|2.2p but such 
that £ := Yl n &N^-( n ) < 1' i- e -' ^(0 i s a sub-probability on N. Then, one may define 
K(n) := K(n)/Ti which is obviously a probability It is easy to realize from Eq. (12, 8ft that 
the Gibbs measure (and free energy) of the model defined starting from K(.) is the same 
as that obtained starting from K(.), provided that h is replaced by h + logS. The case 
where r are the zeros of the SRW on Z d with d > 3 can then be included in our discussion: 
Eq. (|2.2p holds with a = d/2 — 1 and L(.) asymptotically constant. In the following we 
will therefore always assume, without loss of generality, that t is recurrent. 

We conclude this section by listing a couple of examples of (bio)-physical situations 
where disordered pinning models are relevant: 

• Wetting of (1 + 1) -dimensional disordered substrates |17] |21j . Consider a two- 
dimensional system at a first order phase transition, e.g., the 2d-Ising model at 
zero magnetic field and T < T c , or a liquid-gas system on the coexistence line. 
Assume that the system is enclosed in a square box with boundary conditions im- 
posing one of the two phases along the bottom side of the box and the other phase 
along the other three sides. For instance, for the Ising model one can impose + 
boundary conditions (b.c.) along the bottom side and — b.c. along the other ones; 
for the liquid-gas model, one imposes that the bottom of the box is in contact with 
liquid and that side and top walls are in contact with gas. Then, there is necessar- 
ily an interface joining the two bottom corners of the box and separating the two 
phases. At very low temperature, it is customary to describe this interface as a 
one-dimensional symmetric random walk (not necessarily the SRW) conditioned to 
be non-negative, the non-negativity constraint reflecting the fact that the interface 
cannot exit the box. The directed character of the random walk implies in partic- 
ular that one is neglecting the occurrence of bubbles or overhangs in the interface. 
An interesting situation occurs when the bottom wall is "dirty" and at each point 
has a random interaction with the interface: at some points the wall prefers to 
be in contact with the gas (or — phase), and therefore tries to pin the interface, 
while at other points it prefers contact with the liquid (or + phase) and repels the 
interface. Of course, this non-homogeneous interaction is encoded in the charges 
LO n . In this context, the (de)-localization transition is called wetting transition. 
This denomination is clear if we think of the liquid-gas model: the localized phase 
corresponds to an interface which remains at finite distance from the wall (the wall 
is dry), while in the delocalized phase there are few interface- wall contacts and the 
height of the liquid layer on the wall diverges in the thermodynamic limit: the wall 
is wet. It is known that, in great generality |19j . the law of the first return to zero 
of a one-dimensional random walk conditioned to be non-negative is of the form 
(|2,2p with a = 1/2 and L(.) asymptotically constant (this process is transient but 
this fact is not so relevant, in view of the discussion at the beginning of the present 
section) . 

• Formation of loops under thermal excitation and denaturation of DNA molecules 
in the Poland-Scheraga (PS) approximation pjS]. Neglecting its helical structure, 
the DNA molecule is essentially a double strand of complementary units, called 
"bases". Upon heating, the bonds which keep base pairs together can break and 
the two strands can partly or entirely separate (cf. figure below). This separation, 
or denaturation, can be described in the context of our disordered pinning models. 
The set r represents the set of bases whose bond is not broken. In the localized 
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phase t contains O(N) points (N being interpreted here as the total DNA length), 
i.e., corresponds to the phase where the two strands are still essentially tightly 
bound. In the delocalized (or denaturated) phase, on the contrary, only few bases 
pairs are bound. In formulating the PS model, one usually takes a value a ~ 2.12 
(cf. [35J for a justification of this choice) and (in our notations, which are not 
necessarily those of the literature on the PS model) 

L(n) = a for n > 2, 

where a (the cooperativity parameter) is a small number, usually of the order 10~ 5 , 
while L(l) is fixed by the normalization condition XmeN K( n ) = 1- Quenched 
disorder corresponds here to the fact that bases of the different types are placed 
inhomogeneously along the DNA chain. We refer to \22\ Section 1.4] for a very 
clear introduction to the denaturation problem and the Poland-Scheraga model. 
Here we wish to emphasize only that the renewal process r described by such 
a K(.) is not in general the set of returns of a Markov chain, as it happens for 
instance in the case of the wetting model described above. 



Binding energy: Bound base pair 

E n j= (3oj n + h 




Loop Open base pair 



3. The questions we are interested in 

The main questions which will be considered in these notes are the following: 

(1) When is the annealed bound (|2.17j) a good one, i.e., when are quenched and an- 
nealed systems similar? We will see that quenched and annealed free energies never 
coincide, except in the (trivial) case where the annealed free energy is zero (i.e., 
the annealed model is delocalized). However, this does not mean that the solution 
of the annealed system gives no information about the quenched one. For instance 
we will show that, for a < 1/2 and weak enough disorder, the quenched critical 
point coincides with the annealed one. This will be discussed in Section 15.31 

(2) What is the order of the transition? Critical exponents (in particular, the specific 
heat exponent, cf. next section) can be exactly computed for the homogeneous 
model. The Harris criterion predicts that for small (3 critical exponents are those 
of the (3 = (or annealed) model if a < 1/2, and are different if a > 1/2. This is 
the question of disorder relevance, discussed in Sections I5.3ti5.6t 

(3) Truncated correlations functions are known to decay exponentially at large dis- 
tance, in the localized phase. What is the behavior of the correlation length when 
the transition is approached? We will see that, due to the presence of quenched 
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disorder, one can actually define two different correlation lengths. In specific cases, 
we will identify these correlation lengths and give bounds on the critical exponents 
which govern their divergence at h c {(3). 

4. The homogeneous model 

In absence of disorder (/? = 0) the model is under full mathematical control; in particu- 
lar, critical point and the order of the transition can be computed exactly. In this section, 
we collect a number of known results, referring to |22l Chapter 2] for their proofs. 

The basic point is that the free energy F(0, h) is determined as follows [261 Appendix 
A] : if the equation 

e- bn K(n) = e~ h (4.1) 

neN 

has a positive solution b = b(h) > then F(0,h) = b(h). Otherwise, F(0,h) = 0. From 
this (recall the normalization condition Y2neN-^-( n ) = -0> one finds immediately that 
h c (0) = 0. The behavior of the free energy in the neighborhood of h c (0) can be also 
obtained from (|4.ip . Care has to be taken since a naive expansion of left- and right-hand 
sides of (14.11) for b and h small does not work in general. However, this analysis can be 
performed without much difficulty and one can prove the following: 

Theorem 4.1. [221 Th. 2.1] 

(1) If a = 0, F(0, h) vanishes faster than any power of h for h \, 0. 

(2) If < a < 1 then for h > 

F(0,h) = /i 1/Q L(l//i), (4.2) 
where L(.) is the slowly varying function 

L(l/h) = (j^aj) h ~ 1/aR ^ h ) ( 43 ) 

and R a {.) is asymptotically equivalent to the inverse of the map x — > x a L(l/x). 

(3) If a = 1 and XmeN n -^( n ) = 00 then F(Q,h) = hL(l/h) for some slowly varying 
function L{.) which vanishes at infinity. 

(4) If X^neN nK{n) < oo (in particular, if a > 1) 

f(0 ; a) "°e^; <4 - 4) 

In particular, note that in the situation (4), i.e., if r is positively recurrent under P, the 
transition is of first order: the free energy is not differentiable at h c (0) = 0, i.e., the average 
contact fraction has a finite jump in the thermodynamic limit. This is analogous to what 
happens for the Ising model in dimension d > 2: if T < T c and one varies the magnetic 
field H from CP to + , the spontaneous magnetization has a positive jump and the free 
energy is not differentiable. The transition is, on the other hand, continuous (at ieast of 
second order) if P is the law of a null-recurrent renewal r and it becomes smoother as a 
decreases. In thermodynamical language, one can say that the derealization transition 
is of k th order (F(j3, .) is of class C fc_1 but not of class C k ) for a 6 (1/k, 1/k — 1) and of 
infinite order for a = 00 



1 In order to decide between k th and (k + l) th order for a — 1/k one needs to look also at the slowly 
varying function L(.), as is already clear from points (3) and (4) in the case of k = 1. In any case, the 
precise statement is that of Theorem 14.11 
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In the physics literature one introduces usually the specific heat critical exponent v ag^ 

log F(0,h) 

v = 2- hm - — — (4.5) 

h\h c {i3) log(/t - h c {(5)) 

(provided the limit exists) and of course v can depend on (3. From Theorem 14.11 we see 
that, in absence of disorder, 

i/(/? = 0) = 2-max(l,l/a). (4.6) 

In particular, note that v{(3 = 0) > as soon as a > 1/2 (this observation will become 
interesting in the light of the results of Section I5.6[) . 



5. Relevance or irrelevance of disorder? 

We have just seen that the phase transition of the homogeneous pinning model can be 
of any given order - from first to infinite - depending on the choice of K{.) in (|2.2p and, 
in particular, on the value of a. In this section we discuss the effect of disorder on the 
transition and we are primarily interested in the question of disorder relevance. There are 
actually two distinct (but inter-related) aspects in this question: 

Ql does an arbitrarily small quantity of disorder change the critical exponent v (i.e., 
the order of the transition)? 

Q2 does the quenched critical point differ from the annealed one for very weak disor- 
der? 

One expects the answer to both questions to be "no" if a < 1/2 and "yes" if a > 1/2, 
while the case a = a c = 1/2 is more subtle and not clear even heuristically [TTJ [21] (see, 
however, Theorem I5.5H . 

The plan is the following: we will first of all (Section 15. ip make a non-rigorous compu- 
tation, in the spirit of the Harris approach [33], which shows why the watershed value for 
a, distinguishing between relevance and irrelevance, is expected to be a c = 1/2, i.e., the 
value for which the critical exponent v vanishes for the homogeneous model (cf. (|4.6p ). 
Next, in Section [5.21 we prove an upper bound for the free energy which strictly improves 
the annealed bound (|2.17p . In the proof of this bound we introduce the technique of inter- 
polation, by now classical in spin glass theory but sort of new in this context. We would 
like to emphasize that interpolation (and replica coupling, cf . Section 15. 5p techniques have 
proven recently to be extremely powerful in the analysis of mean field spin glass models, 
cf. for instance [32], [I], [39], while their relevance in the domain of disordered pinning 
model had not been realized clearly so far. 

As a byproduct, our new upper bound partially justifies the heuristic expansion of 
Section 15.11 The question of relevance is taken up more seriously in Sections 15.31 to 15.61 In 
the former we will see, among other results, that answers to both Ql and Q2 are actually 
"no" for a < a c . In the latter, on the other hand, we show that critical exponents are 
modified by disorder for a > a c : in particular, we will see that v < whenever (3 > 0. 

In the whole of Section [5] we consider only the case of Gaussian disorder. This allows 
for technically simpler proofs, but results can be generalized for instance to the bounded 
disorder case. 



the symbol v for the specific heat exponent is not standard in the literature, but we have already 
used the letter a for another purpose. The same remark applies to the symbols we use for other critical 
exponents. 
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5.1. Harris criterion and the emergence of a c = 1/2. Let us note for clarity that, 
putting together the discussion of Section 0] and Eq. (12.180 . in the Gaussian case the 
annealed critical point equals h^.{(3) = — /3 2 /2. The first step of our heuristic argument is 
rigorous and, actually, an immediate identity: 

F N (P,h) = F a N {(3,h) + lElog/e^i^"^ 2 / 2 )' 5 ") (5.1) 

where Otv/i := ^^ (.) is just the Boltzmann average for the homogeneous system (cf. 
Eq. (|2.8p ). Identity (|5.ip can be rewritten in a more suggestive way if we recall the last 
equality in (pTTTj) and we let h = h^{/3) + A with A > 0: 

F N (P, h a M + A) = F N (0, A) + R NA (J3) := F N (0, A) + -J-E log ( e ^(K-^/2)^\ (5 . 2) 

iV \ / N,A 

Irrelevance of disorder amounts to the fact that, for (3 sufficiently small, the "error term" 
Rn,a(P) 1S negligible with respect to the "main term" F/v(0,A). As we will see, the 
question is subtle since we are interested in both A and (3 small, and the two limits do not 
in general commute. For the moment, let us proceed without worrying about rigor and 
let us expand naively Rn,a(P) for (3 small and A, N fixed: 

N o2 N 

e E^ l( K-^/2)^^ = ! + £ (/3Wn _ 2 /2) {6n)NA + P_ J- u; n Lo m (5 n 8 m ) NA + 0(/3 3 )(5.3) 

n=l n,m=l 

Expanding the logarithm and using the fact that Ew„ = and E(u; n u; m ) = lr n=m \ one 
has, always formally 

Rna(P) = -|y E 0^,a) + 0(P 3 ). (5.4) 

n=l 

In the limit iV — ► cxd one has by definition of the homogeneous model 

1 N 

lim <^)^ A = lim - J] (6 n ) N A = d A F(0, A). 

n=l 

Since {S n ) N A should not depend on n as soon as 1 <C n <C N, one can expect (actually, 
this can be proven without much difficulty) that 

lim (5 n ) NA = d A F(0,A). (5.5) 

n/A^me(0,l) 

In conclusion, we find 

F((3, h a M + A) = F(0, A) - ^(d A F(0, A)) 2 + 0(/3 3 ). (5.6) 

Even without trying (for the moment) to justify this expansion or to look more closely 
at the A-dependence of the error term 0(/3 3 ), we can extract something important from 
Eq. (ISTB]) . We know from Theorem O that, for a < 1 and A > small, F(0, A) ~ A 1 /" 
which implies (cf. the proof of Eq. (|5.19p for details) that d A F(0, A) ~ A 1 /" -1 . Then we 
see immediately that, indeed, for a < 1/2 

^(d A F(0,A)) 2 «F(0,A) (5.7) 

if A and (3 are small. In terms of the Harris criterion, disordered is said to be irrelevant in 
this case and one can hope that the expansion can be actually carried on at higher orders. 
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For 1/2 < a < 1, however, this is false: even if (3 is small, choosing A sufficiently close to 
zero the left-hand side of (|5.7f) is much larger than the right-hand side. This means that 
"disorder is relevant" and the small-disorder expansion breaks down immediately. The 
same holds for a > 1, when -F(0, A) ~ A and d^Ffo, A) ~ const. The threshold value 
a c = 1/2 is clearly a "marginal case" where relevance or irrelevance of disorder cannot be 
decided (even on heuristic grounds) by a naive expansion in (3. 

The rest of this section will be devoted to give rigorous bases to this suggestive picture. 
As a byproduct we will learn something interesting for the case 1/2 < a < 1: while disorder 
is relevant and changes the exponent v, it modifies the transition only "very close" to the 
critical point (cf. Theorem 15. 3ft . 

5.2. A rigorous approach: interpolation and an improvement upon annealing. 
In Section [2.51 we saw that a simple application of Jensen's inequality implies F(fl,h) < 
F a (f3, h). Here we wish to show that this inequality is strict as soon as disorder is present 
{(3 > 0) and the annealed system is localized. Moreover, we will partly justify the small-/? 
expansion of Section 15.11 for a < 1/2, showing that it provides an upper bound for the 
quenched free energy. 
More precisely: 

Theorem 5.1. @3l Th. 2.6] For every (3 > 0, a > and A > 

F((3, h a M + A) < inf (^f- + F(0, A - (3 2 q)) < F(0, A) = F a ((3, h). (5.8) 
In particular, if < a < 1/2 there exist constants (3q > 0, Aq > such that 

F((3, h a M + A) < F(0, A) - (d A F(0, A)) 2 (1 + 0{f3 2 )) (5.9) 

for (3 < (3q,A < Aq, where 0{(3 2 ) is does not depend on A. On the other hand, if f3 = 
or A < 0, then F{j3, h a c {(3) + A) = F a {(3, h a c {l3) + A). 

About the possibility of pushing the upper bound (|5.9|) to order higher than f3 2 see 
Remark 3.1 in [33]. It is obvious that (|5.9p cannot hold for a > 1/2 since, as already 
observed after Eq. (15. 7p . the right-hand side is negative for A sufficiently small. 

Readers familiar with mean field spin glass models will remark a certain similarity 
between the variational bound (15. 8j) and the "replica symmetric" variational bound [30j 
for the free energy of the Sherrington-Kirkpatrick model. However, we do not see a natural 
way to generalize (15.8P to include "replica symmetry breaking" in analogy with |29j pp. 

Proof of Theorem \5.1\ The proof is rather instructive because it allows us to introduce 
the technique of "interpolation", which will play a major role in the next subsection. We 
start from identity (|5.2p and, for A > 0, q G R and < t < 1, we define 

Rn a(*, A q) :=-E log L^n=^tu> n -tP/2+p\(t-i)]6 n \ (5 10) 

N \ I A,N 

In spin glass language, this would be called an "interpolating free energy" , since by varying 
the parameter t it relates in a smooth way the quantity we wish to estimate at t = 1, 

RnaV = i>M = Rna(0) ( 5 - n ) 

to something easy at t = 0: 

RNA(t = 0, A q) = F N (0, A - P 2 q) - F N (0, A). (5.12) 
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A priori, there is no reason why Rn,a(^j Pi ?) should be any easier to compute for < t < 1 
than for t = 1. What helps us is that the t-derivative of RnA&i P> O) can be bounded 
above by throwing away a (complicated) term which, luckily, has a negative sign. To see 
this we need first of all manageable notations and we will set 



^( r ) e E^=i^vW-t/3 2 /2+/3 2 g(t~l)]<5„^ 



{9{t))n ^ '- " /^^-^/a+^t-i^N ' ~ (5 ' 13) 

\ / AT, A 

for every measurable function g(r). We find then 

d 8 2 ( 1 \ N a N 

^ ' m=l v m=l 

The last term of (|5.14p can be rewritten using the Gaussian integration by parts formula 

E(uf(u))=Ef'(u), (5.15) 

which holds (if u is a Gaussian random variable A/*(0, 1)) for every differentiable function 
/(.) such that lim| a ,|^ 00 exp(— x 2 /2)f(x) = 0. In our case, the function / is of course 
(8m) n A t an d one finds 

3 N P 2 N ( 2 \ 

v m=l m=l ^ ' 

The positive term comes from the differentiation of the numerator of {S m ) N a , (recall the 
definition (|5. 13|) ) and the negative one from the denominator, and we used the obvious 
8m = {&m) 2 - Putting together Eqs. (|5.14p and (|5.16p one has therefore 

j ol 2 Q2 n f . N 2"| fl2 2 

^RNA(t,M = ^-^Y, E {(( S n)NA,t-<l) <^f- ( 5 - 17 ) 

n=l ^ ^ 

At this point we are done: we integrate on t between and 1 inequality (15.17p . we recall 
the boundary conditions (|5.12p and (|5.1ip and we get 

Rna(P) < F n(0, A - 3 2 q) ~ F N (0, A) + (5.18) 

Together with Eq. (|5.ip . taking N ^ oo limit and minimizing over q proves (|5.8p . Let 
us remark that minimizing over q E R or on < q < A//? 2 is clearly equivalent. The 
strict inequality in (|5.8p is just due to the fact that the derivative with respect to q of the 
quantity to be minimized, computed at q = 0, is negative. 

The expansion (|5.9p is just a consequence of (|5.8p . Remark first of all that, at the lowest 
order in 8, the minimizer in (|5.8p is g = (/a := 9aF(0, A). Then, from identity (|4.ip one 
finds that there exist slowly varying functions L^\.),i = 1,2 such that for a < 1/2 and 
A > 

d A F(0, A) = A( 1 - Q )/ Q L( 1 )(1/A), d|F(0, A) = A^^L^l/A). (5.19) 

Let us show for instance the first equality. Differentiating both sides of (|4.ip with respect 
to A one finds 
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Using Theorems IA.1I and I A. 21 one has then, for A — ► (i.e., for F(0, A) — > 0) 

f) j?(n A\ AX0 r(2-a)L(l/F(0,A)) 

9aF( °' A) ~ (l-a)F(0,A)i-« (5 ' 21) 

which, together with (|4.2p . proves the first equality in (|5.19p for a suitable L^'(.). Note, 
by the way, that thanks to (15.191) one has q A < A//3 2 for A, sufficiently small (and 
a < 1/2, of course). Another consequence of (I5TT9|) is that d%F(0, A) is bounded above 
by a finite constant C for, say, A < 1. Then, a Taylor expansion gives 

F(Q, A - 2 q A ) < F(0, A) - 2 (d A F(O, A)) 2 + C0 4 (d A F(O, A)) 2 , 



whence Eq. 

Finally, the statement for = or A < is trivial: for = there is no disor- 
der to distinguish between quenched an annealed free energies, and for A < one has 
F a (0, h a c {0) + A) = which, together with (ETLTD and F((3, h) > 0, implies the statement. 

□ 

5.3. Irrelevance of disorder for a < 1/2 via replica coupling. We want to say first 
of all that, if < a < 1/2 and is sufficiently small (i.e., if disorder is sufficiently weak), 
then h c {p) = h a c {0). Recalling that F a ((3, h a c {(3) + A) = F(0, A), this follows immediately 
from 

Theorem 5.2. [HIE] Assume that either < a < 1/2 or that 

a = 1/2 and ^ n~ l L{n)- 2 < oo. (5.22) 

neN 

Then, for every e > there exist 0o(e) > and Ao(e) > such that, for every < 0o(e) 
and < A < Ao(e), one has 

(1 - e)F(0, A) < F(0, h a c (0) + A) < F(0, A). (5.23) 

Observe that this implies in particular that, under the assumptions of the theorem, the 
exponent v equals 2 — 1/a as in the homogeneous case. Indeed note that, for A small, 

log(l-e) + logF(0,A) > logF(0,h c (0) + A) > logF(0,A) 
log A ~ log A log A 

and the statement follows taking the limit A — ► from definition (|4.5f) of the specific heat 
exponent. 

We will see in Section 15.61 that the same cannot hold for a > 1/2: in that case, v is 
necessarily non-positive in for the quenched system presence of disorder, while it is positive 
for the annealed system. One could therefore think that quenched and annealed behaviors 
are completely different. This is however not completely true. Indeed, the next theorem 
shows that F(0,h) and F a (0,h) are very close, provided that 1/2 < a < 1 if one is not 
too close to the critical point. More precisely one has 

Theorem 5.3. Assume that 1/2 < a < 1. There exists a slowly varying function L{.) 
and, for every e > 0, constants ai(e) < oo and Ao(e) > such that, if 

a 1 (e)0 2a ^ 2a - 1 ^L(l/0) < A < A (e), (5.25) 

the inequalities (15.23)) hold. 
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To see more clearly what this says on the relation between quenched and annealed 
critical points, forget about the slowly varying functions; then, Theorem 15.31 implies 

0<h c (P)-h a c (f3)<P 2a/(2a ~ 1} - 
Since 2a/ (2a — 1) > 2, this shows in particular that 

lim M§ = 1. (5.26) 

Remark 5.4. Theorem 15.31 was proven in [U Th. 3] and then in [43|, Th. 2.2]. The two 
results differ only in the form of the slowly varying function £(.). In general, the function 
L(.) which pops out from the proof in |43j Th. 2.2] is larger (i.e., worse) than that of [U 
Th. 3]. 

Finally, we consider the "marginal case" a = a c = 1/2 and X^n(-^( n ))~ 2n_1 = 00 • This 
is the case, for instance, if P is the law of the returns of a one-dimensional symmetric 
random walk, where L(.) is asymptotically constant, as mentioned in Section [2.61 As we 
mentioned, this case is still debated even in the physical literature. The "most likely" 
scenario [T7j is that disorder is "marginally relevant" in this case: h c (0) ^ h^(/3) for 
every positive j3, but the two critical points are equal at every order in a weak-disorder 
perturbation theory. Other works, e.g. |21j . claim on the other hand that disorder is 
irrelevant in this situation. 

What one can prove for the moment is the following: 

Theorem 5.5. [H H3] Assume that a = 1/2 and X^neN n " 

_1 L(n)- 2 = oo. Let £(.) be the 

slowly varying function (diverging at infinity) defined by 

N 1 

n=l v ; 

For every e > there exist constants 02(e) < 00 and Ao(e) > such that, ifO < A < Ao(e) 
and if the condition 

1 f ^f a 2 (e)\logF(0,A)\ \ 

p*<»M{ FftA) J (5 ' 28) 

is verified, then Eq. (|5.23|) holds. 

Remark 5.6. To be precise, in the statement of [3] Th. 4] the condition (|5,28p is replaced 
by a different one (essentially, the factor |log-F(0,A)| in the argument of £(.) does not 
appear). In this sense, the condition (|5.28p under which we prove here (|5.23p is not the 
best possible one. However, for many "reasonable" and physically interesting choices of 
L(.) in (|2.2p . Theorem 15.51 and Theorem 4 of [3] are equivalent. In particular, if P is 
the law of the returns to zero of the simple random walk {S n } n >o in one dimension, i.e. 
r = {n > : $2™ = 0}, in which case L(.) and L(.) are asymptotically constant and 
£(N) ~ 03 log N, one sees easily that (|5.28|) is verified as soon as 

A>a 4 (e)e~^ 1 , (5.29) 

which is the same condition given in [3]. 

Note, by the way, that in this case the difference h c ((5) — h^((3) vanishes faster than any 
power of j3, for (3 \ 0. This confirms the fact that, even if the two critical points can be 
different, they cannot be distinguished perturbatively. 
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5.4. Some open problems. The results of previous section, while giving rigorous bases 
to predictions based on the Harris criterion, leave various intriguing gaps in our compre- 
hension of the matter. Let us list a few of them, in random order: 

• Let a < 1/2. Does there exist a j3 c < oo such that h c (@) ^ h^(P) for j3 > /3 C ? If 
yes, how smooth is h c ((3) at /3 C ? Does v equal 2 — 1/a also for j3 large? 

• Again, let a < 1/2 and look at Eq. (|5.9p . Is it true that 

F(P, + A) > F(0, A) - ^(d A F(0, A)) 2 (l + 0(/? 2 ))? 

• Under the assumptions of Theorems 15.31 or 15. b\ does there exist positive values of 
(3 for which quenched and annealed critical points coincide? It is sort of reasonable 
to conjecture that the answer is "no", at least for a > 1/2. 

The reader might be tempted to think that such questions should be easy to answer 
numerically. If so, he should have a look at Ref. [12] where one gets an idea (in the 
context of random heteropolymers at selective interfaces) of why numerical tests become 
extremely hard in the neighborhood of the critical curve. 

Remark 5.7. Between the time these notes were written and the time they were pub- 
lished, the above open problems have been to a large extent solved. In particular: 

• in Ref. |44J it was proven that for every a > 0, if (3 is large enough and, say, uj is 
Gaussian, then h c {(3) ^ h^(f3). 

• The question posed in open problem (2) has been answered positively in Ref. |28j . 
although in a slightly weaker sense. 

• In Ref. [16] it was proven that as soon as a > 1/2 and (3 > one has h c (j3) ^ h^(f3). 

5.5. Proof of Theorems I5.2ti5.5l We follow the approach of [43J which, with respect 
to that of [1], has the advantage of technical simplicity and of being closely related to the 
interpolation ideas of Section 15.21 On the other hand, we encourage the reader to look 
also at the methods developed in [3], which have the bonus of extending in a natural way 
beyond the Gaussian case and of giving in some cases sharper results (cf. Remarks 15.41 
and 15.61 above) . 

A natural idea to show that quenched and annealed systems have (approximately) the 
same free energy is to apply the second moment method: one computes E(£/jv(/3, h)) and 
E((Z/v(/?, h)) 2 ) and if it happens that the ratio 

\EZ N ,UP,h)] 2 , 

ekzjvGs./i)) 2 ] 1 ' ' 

remains positive for N —* oo, or at least it vanishes slower than exponentially, it is not 
difficult to deduce that F(/3, h) = F a (f3, h). This approach has turned out to be very pow- 
erful for instance in controlling the high-temperature phase of the Sherrington-Kirkpatrick 
mean field model in absence of magnetic field [40| Ch. 2.2]. However, this simple idea does 
not work in our case and the ratio (|5.30p vanishes exponentially for every j3, A > 0. This 
is not surprising after all, since we already know from Theorem 15.11 that quenched and an- 
nealed free energy do not coincide. There are two possible ways out of this problem. One 
is to perform the second moment method not on the system of size N but on a smaller sys- 
tem whose size N(A) remains finite as long as A is positive and fixed, and diverges only for 
A — > 0. If N(A) is chosen to be the correlation length of the annealed system, one can see 
that on this scale the ratio (|5.30p stays positive, so that F N ^(/3, h^((3)+A) ~ F N ^(0, A). 
One is then left with the delicate problem of glueing together many blocks of size N(A) 
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to obtain an estimate of the type F(@, h^{(5) + A) > (1 — e)F(0, A) for the full free energy. 
This is, in very rough words, the approach of Ref. [4J. The other possibility, which we 
are going to present, is to abandon the second moment idea in favor of a generalization of 
the replica coupling method [31J [43J. This method was introduced in [31] in the context 
of mean field spin glasses and gives a very efficient control of the Sherrington-Kirkpatrick 
model at high temperature {(5 small), i.e., for weak disorder, which is the same situation 
we are after here. 

The two methods are in reality not orthogonal: they share the idea that the important 
object to look at is the intersection of two independent renewals To see why 

this quantity arises naturally, let us compute the second moment of the partition func- 
tion. If tM,t( 2 ) are independent renewal processes with product law P® 2 (.), recalling the 
definition A = h + /3 2 /2, one can write 

E((Z N ,UP,h)f) = EE 02 (e E "- ( ^ +h){1 { «-< 1 )> +1 ( -^ 2 »} ) l {iVer(1)} l {JVeT(2)} ) (5.31) 



E® 2 



3 A(|r( 1 )n{i,...,7V}|+|r( 2 )n{i,...,iV}|)+/3 2 |r( 1 )nr( 2 )n{i,...,iV}| 1 , 



Considering also that 

12 _ P ,®2 / A|(r( 1 )n{l,...,JV}|+|r( 2 )n{l,...,7V}|)- 



tW 7 (R h\\ z — TP®^ / t\\ (r^ T i|l,...,iv } + r v T l| l,...,iv } J -i -. 

one sees that the ratio (I5.30P depends on the typical number of points that and 
have in common up to time N. One sees also why this ratio has to vanish exponentially 
N — ► oo: as long as A > the renewals r^, with law modified by the factor exp(A|r( 1 ) n 
{1, . . . , iV}|), are finitely recurrent and therefore will have a number of intersections in 
{1, . . . , N} which grows proportionally to N. 

Proof of Theorem \5.2{ The second inequality in (15.231) is just Eq. (12.170 . As for the 
first one, let A > and recall identity (j5.2f) . Define, in analogy with (I5.10|) . 

R NA (t,/3) := lllog/ e ^»=x(^«-^ 2 /2)5n\ (5 .3 2) 
N \ I A,N 

for < t < 1 (to the purpose of Theorem 15.21 we do not need the variational parameter q) 
where the measure {-) NA was defined after Eq. (|5.ip . Observe that 

Rn,a(°,P)=° ( 5 - 33 ) 

while 

RnaO->0)=Rn,a(P)- ( 5 - 34 ) 
As for the t-derivative one finds (just take (|5.17p and put q = 0): 

' /^ me E^=i(/V^n-t/3 2 /2)<C 



m=l 



A,N 



3 En=i(/Vfc^-*/3 2 /2)5 n 



(5.35) 



A,N 



Recall definition (|5.13p (specialized to the case q = 0) of the random measure (.) N A t and 
let (-}% 2 ^ t be the product measure acting on the pair (t^\t^), while 5$ '■= Ijnert*)}- 
Note that the two replicas T®,i = 1,2 are subject to the same realization uj of disorder. 
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Then, one can rewrite 

?2 N 



m=l 

Since we need a lower bound for Rn A (P) to prove the first inequality in (|5.23p . the fact 
that this derivative is non-positive seems to go in the wrong direction. Let us not lose 
faith and let us define, for A > 0, 

2N & \ /n,A 

_L E log / e E^=iC3^n-^ 2 /2)(^ l3 +^ 2) )+A^ ELi S ( n ] S^ 



(5.37) 



2N " \ I N,A 

where the product measure (.)® 2 a ac ts on the pair 

(rW,r (2) )- The index "( 2 )" refers to 
the fact that this quantity involves two copies (replicas) of the system. Observe that we 
are letting the two replicas interact through a term which is positive, extensive (i.e., of 
order N) and closely related to what appears in the right-hand side of Eq. (|5.36p . Note 
also that 



< } A (0,A,/3) 2V 

while the factor 2 in the denominator guarantees that 

B$ A (t,0,P)=R NA (t,P). (5.39) 

Again via integration by parts (the computation is conceptually as easy as the one which 
ledtoEq. ([537]) ). 



1 " 1 A/3 2 E£ =1 *iM 2) \® 2 



log ( e 



AT, A 



(5.38) 



dt W < A 



2iV ^ 



m=l 



(5.40) 



AT, A 



4N 



m=l 



/((^ ) +4 2) )e^^ (1) - (2) ) 



V 



,.ffjv(f,A ) y8;T( :l ),r( 2 )) 



AT, A 




. 2 iv (*<MW^V<*>)\® 2 

P m \ / N,A 



- 2iV ^ 

m=l 



,H N {t,\,P;TW,TW) 



— i?^ A (t, A,/3). 



AT, A 



This can be rewritten as 



dt 

which implies that, for every < t < 1 and A 

? ( 2 ) u \ a\s d(2) 



R%> A (t,\-t,0) <0 



R%> A (t,\,P) < R% A (0,\ + t,p). (5.41) 

Going back to Eqs. (|5.35|) and the last equality in (|5.4(jp and using the fact that for every 
convex function ip(.) one has xip'(0) < ip(x) — ip(0) one finds 



-^(-R NA (t,p)) = A R%] A (t,\,P) 



RP A (t,2-t,P)-R% A (t,0,P) 



(2) 



< 



A=0 



2-t 



(5.42) 
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Finally, using monotonicity of R N \(t, A, ($) with respect to A and (|5.39p . one obtains the 
bound 

(-R NA (t,/3)) < RP A (O,2,p) + (-R NA (t,0)), (5.43) 

where we used (|5.4ip and the fact that 2 — t > 1 (of course, we could have chosen 1 + rj — t 
instead of 2 — t for some rj > in (|5.42p and the estimates would be modified in a 
straightforward way). We can now integrate with respect to t between and 1 this 
differential inequality (or use GronwalFs Lemma, if you prefer) and, recalling Eqs. (|5,34l) 
and (|5.33p . we obtain 

- (e - 1)2$^ (0,2,0) < R NA (P) < 0. (5.44) 

Before we proceed, we would like to summarize what we did so far. To prove Theorem 
I5.2l we need the lower bound limjv-»oc Rn A (P) > —eF(0, A) but, as in Section I5T21 it seems 
that the interpolation method gives rather upper bounds on Rn A (P)- Then, through the 
replica coupling trick we transferred this problem into the problem of proving an upper 

(2) 

bound for a quantity, R N A (t, A, /3), which is analogous to Rn A (P), except that it involves 
two interacting copies of the system. Moreover, by throwing away a (complicated, but 
with a definite sign) term in Eq. (I5,40p . we reduced to the problem of bounding from 

(2) 

above R N A (t, A, (3) computed at t = 0. In other words, we replaced the task of estimating 
from below Rn A (P) with that of estimating from above a quantity which involves no 
quenched disorder, and which is therefore easier to analyze. While this procedure might 

(2) 

look a bit magic, the basic underlying idea is the following. R N A (i, A, f3) is obviously 

non-decreasing as a function of A. Suppose however that, for some A > 0, R N A (t, A, (3) 
is not very different from the value it has at A = (of course, proving this amounts to 

(2\ 

proving an upper bound on R N A (t, A, /?).) Then, looking at the definition (I5.37p . this 

means that the cardinality of the intersection n n {1, . . . , N} is typically not large 
and this, through Eqs. (|5.33p . (|5.34p and (|5.36p implies a lower bound on Rn A (@)- 
Let us now restart from (|5.44p and note that 

<>,2,/3) = -F N (0, A) + ^logE^ 2 (e^ESLx^^+AESLx^-H^)^)^)" 



< -F N (0, A) + ^ ( °' gA) + — log E 02 ( e^ 2 ESU &® 
q 2Np V 



(5.45) 



where we used Holder's inequality and the positive numbers p and q (satisfying 1/p+l/q 
1) are to be determined. Taking the thermodynamic limit, 

lim sup R%\{0, 2, 0) < lim sup log E 02 ( e 2 ^ ^n=i ®& 

But we know from the expression (|4.2p of the free energy of the homogeneous system and 
from the property (|2.3p of slow variation that, for every q > 0, 

l im ffigfO = (5.47) 
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Therefore, taking q = q(e) sufficiently" close to (but strictly larger than) 1 and Ao(e) > 
sufficiently small one has, uniformly on (3 > and on < A < Ao(e), 

limsup!2SJ ) A (0,2,/3) < — 6 —F(0, A) + lim sup J , . log E® 2 f e 2 ^ 2 E»=i <5 " )(5 " )N )(5.48) 
at—™ ' e - 1 a/ - — >nn 2Nn(e) \ / 



Of course, p(e) := q(e)/(q(e) — 1) < oo as long as e > 0. Note that, in view of (|5.44j) . 
Theorem 15.21 would be proved if the second term in the right-hand side of (|5.48p were 
zero. Up to now, we have not used yet the assumption that a < 1/2 or that (|5.22p holds, 
but now the right moment has come. The way this assumption enters the game is that it 
guarantees that the renewal t^ 1 ' n is transient under the law P® 2 . Indeed, 

E^ 2 I £ l ner(1)nT ( 2 ) I = £ P(n G r) 2 < oc (5.49) 



in>l / n>l 

since, as proven in [18] . 

P("^)"~°°77T^- gS "' ( ' r ° ) f , !,-, ("0) 
L{n)n L a 7r L{n)n L a 

Actually, Eq. (|5.50p holds more generally for < a < 1 and we will need it to prove 
Theorems 15.31 and 15.51 

Transience and renewal properties of the process of n implies that 

p®2(| r (l) n r (2)| > jfe) < (i _ c )* (5.51) 

for some < c < 1: after each "renewal epoch", i.e., each point of rW n t^ 2 \ the 
intersection renewal has a positive probability c of jumping to infinity. Therefore, there 
exists 0i > such that 



N 

for every (3 2 p(e) < (3f. Together with (|5.48|) and (|5.2|) . this implies 

F(P,h a c (P) + A) > (l-e)F(0,A) (5.53) 

as soon as 1 < (3 2 (e) := (3 2 /p(e), and therefore the validity of Theorem 15.21 □ 

Proof of Theorem 15.31 In what follows we assume that A is sufficiently small so that 
F(0, A) <C 1. For simplicity of exposition, we assume also that L(.) tends to a positive 
constant L(oo) at infinity (for the general case, which is not significantly more difficult, 

cf. [13]). 

If we try to repeat the proof of Theorem 15.31 in this case, what goes wrong is that 
the intersection D is now recurrent, so that (|5,52p does not hold any more. The 
natural idea is then not to let -/V tend to infinity at A fixed, but rather to work on a 
system of size iV(A), which diverges only when A — * 0, i.e., when the annealed critical 
point is approached. In particular, we let iV = N(A) := c| log F(Q, A)|/F(0, A) with c > 
large to be fixed later. Note also that this choice of N(A) is quite similar to that made 
in [3], where one applies the second moment method on a system of size c/F(Q, A) with 
c large. This choice has a clear physical meaning: indeed, we will see in Section [6] that 
the correlation functions of the annealed system decay exponentially on distances of order 
1/.F(0, A) (the logarithmic factor in our definition of N(A) should be seen as a technical 
necessity). 



supE 552 (eM^E^s^ < x {5 52) 
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By the superadditivity property (|2.11|) we have, in analogy with (|5.1|) . 

F(/3, -P 2 /2 + A) > F N(A) (0, A) + R N(A)A (f3). (5.54) 

To prove Theorem 15.31 we need to show that the first term in the right-hand side of 
(|5.54p is essentially F(0, A), while the second is not smaller than —eF(0, A), in the range 
of parameters determined by condition (|5.25p . The first fact is easy: as follows from 
Proposition 2.7 of [25], there exists S (0, oo) (depending only on the law K(.) of the 
renewal) such that 

loe N 

F N (0,A)>F(0,A)-a 6 -^- (5.55) 

for every N. Choosing c = c(e) large enough, Eq. (15.550 implies that 

F MA) (0,A)>(l-e)F(0,A). (5.56) 
As for RmM,A(P), w e have from Eqs. (|5.44|) and (|5.45|) 



1 .logE^^E^if'^^ 



(5.57) 



2N(A)p 

where we used Eqs. (|5.56p and (|2.11|) to bound — (l/q)F N ^ (0, qA) + i*V(A)(0, A) from 
below. Choosing again q = q(e) we obtain, for A < Ao(e), 

R »™*W > —2eF(0, A) - logE® 2 (e*^E^tfW . (5 . 58) 



e-1 - v ' ' 2N(A)p{e) 

It was proven in [H Lemma 3] and [43\ Section 3.1] that if 1/2 < a < 1 there exists 
a-j =€ (0, oo), which depends in particular on L(oo), such that for every integers N and k 

P m (| r( i) n r(2) n{i JV>|>*)<(i- , (5.59) 

which should be compared with (|5.5ip . valid for a < 1/2. Thanks to the geometric bound 
(|5.59p we have 

/N(A) \ 

E^M^E^WS ) = ^ P ®2 lj- S m S (2) =k U^* (5.60) 

fc>0 \ 71=1 / 



whenever 



1 ' AT(A) 2 "- 1 



e 2 /? 2 ^) ( i a S ] < i 



and this is of course the case if 



=^ 2 p(e) ( 1 a _l ) < ( J ) ( 5 61 ) 

V iv(A) 2 "- 1 ) - v 2at(a) 2q - 1 y • 1 J 



At this point, using the definition of N(A) and point (2) of Theorem 14. 1} it is not difficult 
to see that there exists a positive constant as(e) such that (|5.6ip holds if 

A (2a-l)/a 

72a-l - 



f3 2 p(e) < a 8 (e) - — —5— T - (5.62) 



|logF(0, A) 
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Condition (|5.62|) is equivalent to the first inequality in (|5.25|) , for a suitable choice of a\ (e) 
and L(.). As a consequence, for N(A) sufficiently large (i.e., for A sufficiently small) 

2iV(AM e ) i0gJij \ e J -2 C ( e )p( e )|logF(0,A)| 1Og V a 7 

Recalling Eq. (|4.2p one sees that, if c(e) is chosen large enough, 

1 dogE^ 2 (eMWXZLVWW) < eF(0, A). (5.64) 



2N{A)p(e) 

Together with Eqs. (|5.54p . (|5.56p and (|5.58p . this concludes the proof of the theorem. □ 

Proof of Theorem 15.51 The proof is almost identical to that of Theorem 15.31 and up 
to Eq. (|5,58p no changes are needed. One has however to be careful with the geometric 
bound (|5.59p : in this case, it is not sufficient to replace a by 1/2, since the behavior at 
infinity of the slowly varying function L(.) in (|2.2p is here essential. The correct bound in 
this case is (cf. [U Lemma 3] and [131 Sec. 3.1]) 

N \ / \ k 



P 



i2 



for every N, for some ag > 0. We recall that £(.) is the slowly varying function, diverging 
at infinity, defined by (|5,27p . In analogy with Eq. (j5.60|) one obtains then 

^MW^iV^e^ < A e 2^(e) A _ ^ (5.66) 

whenever the right-hand side is positive. Choosing 02(e) large enough one sees that if 
condition (|5.28p is fulfilled then 

and, in analogy with (|5.63|) . 

2iV(AM e ) i0gJlj \ 6 J -2c(eM6)|log J F(0,A)| 1 ° g V a 9 

From this estimate, for c(e) sufficiently large one obtains again (|5.64p and as a consequence 
the statement of Theorem 15.51 □ 

5.6. Smoothing effect of disorder (relevance for a > 1/2). Section [5.31 was devoted 
to showing that, for a < a c , (weak) disorder is irrelevant, in that it does not change the 
specific heat exponent v and in that the transition point coincides with the annealed one 
as long as (3 is small. We saw also that for a c < a < 1 quenched and annealed free 
energies and critical points are very close (Theorems 15.31 and 15.51) . This might leave the 
reader with the doubt that disorder might be irrelevant in this situation too. The purpose 
of the present section is to show that this is not the case. 

We start by recalling that via Theorem [4J] and (12.17j) we know that F((3, h%((3) + A) < 
^max(i/a,i) _ This bound is however quite poor: if we go back to (|5.8p and we choose 
q = A/(3 2 we obtain 

F(AW)+A)<^ (5.69) 
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which is better, for A small and a > 1/2. The point is however that, since one expects that 
h%(P) ^ h c {(3) in this situation, (|5.69)) does not say anything about the critical behavior of 
the quenched system; for this, we would need rather an upper bound on F(/3, h c (f3) + A). 
This is just the content of the following result, which we state in the case of Gaussian 
disorder: 

Theorem 5.8. [26l EZ] For every (3 > 0, a > and A > one has 

F(AM/3) + A)<%±^A 2 . (5.70) 



Remark 5.9. Theorem 15.81 actually holds beyond the Gaussian case; for instance, in the 
case of bounded variables u n . In this case the statement has to be modified in that the 
factor 2 in that the denominator in the right-hand side of (|5.70p is replaced by c := c(P), 
a constant which depends only on the disorder distribution P, and the results holds only 
provided A is sufficiently small: A < Ao(P), see |26| . 

Remark 5.10. An obvious implication of Theorem 15.81 is that v < as soon as (3 > 0. In 
this sense, this result is much reminiscent of what was proven in [131 [T4"j about the specific 
heat exponent for the nearest-neighbor disordered Ising ferromagnet. 

In particular, Theorem 15.81 shows that the specific heat exponent is modified by an 
arbitrary amount of disorder if a > a c : the phase transition is smoothed by randomness if 
a > a c and becomes at least of second order (the effect is particularly dramatic for a > 1, 
where the transition is of first order for (3 = 0). 

It is also interesting to compare Theorem 1 5 . 8 1 wit h the celebrated result by M. Aizenman 
and J. Wehr [2] which states that first order phase transition in spin systems with discrete 
spin-flip symmetry are smoothed by disorder as long as the spatial dimension verifies 
d < 2, while the same holds for d < 4 if the symmetry is continuous. 

A less obvious consequence of Theorem 15.81 is the following: 

Theorem 5.11. |41j Let (3 > and < a < oo. There exists c > such that 

lim EP} h J (/3) (Jr n{l,...,Af}| > cAf 2/3 logiV) =0. (5.71) 
Moreover, under the assumptions of Theorem \5.2\ for (3 sufficiently small 

lim EPj h J (/3) (\m{l,...,N}\ >cN 2a ^ 1+a hogN) =0. (5.72) 

This result should be read as follows. The fact that the transition is at least of second 
order in presence of disorder implies already that the Gibbs average of the contact fraction 
defined by (|2.14p tends to zero in the thermodynamic limit at the critical point. The 
additional information provided by Theorem 15.111 are finite- N estimates on the size of 
r n {1, . . . , N} at criticality. Whether the exponent 2/3 in Eq. (|5.7ip is optimal or not is 
an intriguing open question. 

Theorem 15.111 was proven in [41 (together with more refined finite-size estimates on 
E A>d r n {!> • • • i N }\) for h-h c ((3) going to zero with N), apart from Eq. (j57?2|) which 
is a consequence of [41| Th. 3.1] plus Theorem 15.21 (cf. also Remark 3.2 in |41j). 

^ Theorem 3.1 in [IT] is formulated in the case of bounded random variables u n , but it generalizes 
immediately to the Gaussian because the basic ingredient one needs is the concentration inequality [411 
Eq. (5.2)], which holds in the case of Gaussian randomness as well. 
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Proof of Theorem 15,81 (sketch ) For a fully detailed proof we refer to . In the case of 
Gaussian disorder a simpler proof is hinted at in [27] and fully developed in [221 Section 
5.4]. 

Here we give just a sketchy idea of why the transition cannot be of first order when 
/3 > 0. Assume by contradiction that 

F((3, h c (j3) + A) ~ cA for A -> + , (5.73) 

and consider the system at the critical point (f3,h c (/3)). Divide the system of size ./V 
into N/M blocks Bi of size M, with the idea that 1 « M « N. For a given real- 
ization of lo mark the blocks where the empirical average of u, i.e., (1/M) ^] ngB , u> n 
equals approximately A//?. By standard large deviation estimates, there are typically 
Nmarked '■= (N/M)e~ MA ^ 2/3 ) such blocks, the typical distance between two successive 
ones being D typ := Me MA2 /(2/3 2 ). it 

is a standard fact that if we take M IID standard 
Gaussian variables and we condition on their empirical average to be 5, for M large they 
are (roughly speaking) distributed like IID Gaussian variables of variance 1 and average 
5. Therefore, in a marked block the system sees effective thermodynamic parameters 
(P e ff,h e jf) := (f3,h c ({3) + A). Now we want to show that the assumption (|5.73p leads 
to the (obviously false) conclusion that F(j3,h c (j3)) > 0. Indeed, let be the set of r 
configurations such that: 

• there are no points of r in unmarked blocks 

• the boundaries of all marked blocks belong to r. 

Note that depends on disorder through the location and the number of marked blocks, 
and that there is no restriction on r inside marked blocks. One has the obvious bound 

F N (fi,h c {0)) > ^ElogE [e^^ +h ^l {T&Shj] 5 N ) . (5.74) 

But due to the definition of the set , the restricted free energy in the right-hand side of 
(I5.74P gets (for M large) a contribution M mar ked x (M/N)F(/3,h c (j3) + A) from marked 
blocks, and an entropic term J\f m arked/Nx log K(D typ ) from the excursions between marked 
blocks. Summing the two contributions, recalling the asymptotic behavior (|2.2I) of if (.), 
the expression of N mar ked an d F>t yv and taking the 7V — ► oo limit at M large but fixed one 
obtains then 

F((3, h c ((3) + A) > e - MA2 /(^ 2 ) ( F {^ h c {0) + A) - (1 + a)^j . (5.75) 

Since the left-hand side of (|5,75p is zero, for A small and (3 > this inequality is clearly 
in contradiction with the assumption (|5.73p that the transition if of first order (actually, 
even with the assumption F(f3, h c ((3) + A) ~ cA y with y < 2). 

6. Correlation lengths and their critical behavior 

From certain points of view, the localized region C is analogous to the high-temperature 
phase of a spin system. Indeed, in this region one can prove typical high-temperature 
results like the following: free energy fluctuations are Gaussian on the scale 1 / y/N [31 [25] , 
the infinite- volume Gibbs measure is almost-surely unique and ergodic [9], the free energy 
is infinitely differentiable, finite-size corrections to the infinite volume free energy are of 
order 0(1/N), and truncated correlation functions decay exponentially with distance [25J. 
In this section we concentrate on the last point, which turns out to be more subtle than 
expected, in particular when one approaches the critical line. 
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In this section we assume that the random variables oo n are bounded, because the 
results we mention have been proved in the literature under this assumption. They should 
however reasonably extend to more general situations, for instance to the Gaussian case. 

In the following, Pbo^w(-) wu l denote the infinite-volume Gibbs measure, defined as 
follows: first of all we modify definitions (|2.6p and (|2.8j) replacing X^n^i ( A^n + h)S n by 

LiV/2j 

[Pu n + h )5 n , 

n=-[N/2\ 

where {w n } n gz are IID random variables, and then for a local observable /, i.e., a function 
of r which depends only on r D 7 with 7 a finite subset of Z, we let 

E&i,(/) := Hm (/). (6-1) 

Existence of the limit, in the localized phase, for almost every disorder realization is proven 
in (cf. also [9], where a DLR-like point of view is adopted) @ 

The definition of the correlation length £ contains always some degree of arbitrariness, 
but conventional wisdom on universality states that the critical properties of £, close 
to a second-order phase transition, are insensitive to the precise definition. There is 
however a subtlety: in the case of disordered systems there are two possible definitions of 
correlation lengths, which have no reason to have the same critical behavior. Remaining 
for definiteness in the framework of our disordered pinning models, one can first of all 
define a (disorder-dependent) two-point function as 

CUkJ) := P^(k G t\£ G t) - P&JJ* G r). (6.2) 

In words, C UJ (k,£) quantifies how much the occurrence of I G r influences the occurrence 
the event fc G r. It is then natural to define a correlation length £ as 

J :=- lim ^ log|C (A;, 0)|, (6.3) 

provided the limit exists. Note that £ depends on (/?, /i) and, in principle, on u. One can 
however define a different correlation length, £ av , as 

J_:=- lirnilo g E|afc,0)|. (6.4) 

In other words, £ (respectively, £ a,; ) is the length over which the two-point function (re- 
spectively, the averaged two-point function) decays exponentially. For simplicity, we will 
call £ the typical (or quenched) correlation length, and £ av the average correlation length, 
although it is important to keep in mind that £ av is not the disorder-average of £ (indeed, 
in Section I6.3I we will see an example where £ is almost-surely constant but £ 7^ £, av ). 
It is interesting that in the case of the one-dimensional quantum Ising chain with ran- 
dom transverse field studied in [20], the two correlation lengths are believed, on the basis 
of a renormalization group analysis, to diverge at criticality with two different critical 
exponents. 

4 One might give a different definition of the infinite-volume Gibbs measure, considering the original 
system (|2.8p defined in {1, ... , iV} and taking a the N — > 00 limit of the average of local functions of r n J, 
with / a finite subset of N. In other words, with the first procedure, Eq. (|6.f ) , we are looking at the system 
in a window which is situated in the bulk, very far away from both boundaries. On the other hand, the 
second procedure is relevant if one wants to study the system in the vicinity of one of the two boundaries 
(and very far away from the other one). 
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A simple application of Jensen's inequality shows that £ av > £. This inequality can be 
interpreted on the basis of the following intuitive argument. Divide all possible disorder 
realizations into sets A m where the empirical average of to in the region {1, . . . , k} is ap- 
proximately m. Of course, for m ^ A m is a large deviation-like event of probability 
~ exp(— km 2 /2). Conditionally on A m , the system sees a defect line which is more at- 
tractive (if m > 0) or more repulsive (if m < 0) than it should and therefore it is more 
localized (resp. more delocalized) in this region than in the rest of the system. Therefore, 
conditionally on A m , we can expect that C w (/c,0) behaves like exp{—k/^{j3^h + f5m)). In 
other words, we can argue that (looking only at the exponential behavior) 

EC u {k,0)~ J dme-^^e-'^^+W ~ e ^« m {-m 2 /M/$(P+Wl ( 6 . 5 ) 

for k large. Since £ should diverge when the critical point is approached, it is reasonably 
decreasing in h so that the value of m which realizes the maximum is strictly negative. On 
the other hand, when we take the limit without disorder average as in (I6.3p . the events A m 
with m ^ cannot contribute, i.e., almost surely they do not occur for k large enough, as 
follows from the Borel-Cantelli lemma. 

6.1. Correlation length of the homogeneous model. In the homogeneous case, (3 = 
0, the infinite- volume Gibbs measure can be explicitly described (cf. \22\ Th. 2.3]): under 
Poo l (-)j is a homogeneous^, positively recurrent (for h > h c (0) = 0) renewal on Z such 
that 

P^(inf{/fc > : k G r} = n\0 G r) = K{n)e~ F ^ h)n e h =: K h {n) (6.6) 

and 

P°*(n G r) ' 



Note that K^{-) is a probability on N (cf. Eq. (|4.ip and the discussion after it) with an 
exponential tail. What we are interested in is the precise large-n behavior of 

P°*(n G r|0 G r) ' 



E m£ N mK h {m) 

i.e., a refinement of the renewal theorem (which simply states that this quantity tends to 
zero for n — > oo). 

Let us for a moment widen our scope and consider a homogeneous, positively recurrent 
renewal, with law P, such that the law of the distance between two successive points, 
denoted by K (.), has exponential tail: say, 

lim - log Kin) = -z < 0. (6.7) 

We do not require for the moment that K(.) is given by (|6.6|) with K{.) in the class (|2.2|) . 
It is known (cf. for instance [U Chapter VII. 2] and [38]) that, under condition (16. 7p . there 
exist r > and C < oo such that 



P(n G t|0 G t ' 



E meN mK(m) 



< Ce~ rn . (6.8) 



^That is, its law is invariant under translation on Z. For instance, P^ h (n, m G r) = P^/ 1 (n + k,m + k G 
t) for every k £ Z. 
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However, the relation between z and the largest possible r in Eq. (|6.8j) . call it r max , is 
not known in general. A lot of effort has been put by the queuing theory community in 
investigating this point, and in various special cases it has been proven that r max > z (see 
for instance [7J, where power series methods are employed and explicit upper bounds on 
the prefactor C are given). In even more special cases, for instance when P is the law 
of the return times to a particular state of a Markov chain with some stochastic ordering 
properties, the optimal result r max = z is proved (for details, see [SHUT], which are based 
on coupling techniques). However, the equality r max = z cannot be expected in general. 
In particular, if K (.) is a geometric distribution, 

^ „-nc 

K(n) = - 



c - 1 

with c > 0, then one sees easily that the left-hand side of (|6.8p vanishes for every n £ N so 
that r max = oo, while z = c. On the other hand, if for instance K(l) = K(2) = 1/2 and 
K{n) = for n > 3, then z = +oo while r is finite. These and other nice counter-examples 
are discussed in [7J. 

In view of this situation, it is highly non-trivial that, restricting to our original class of 
renewals, the following holds: 

Theorem 6.1. [23J Let Kh{.) be given by (16. 6j) with K(.) satisfying (|2.2[) for some a > 
and slowly varying L(.). Then, there exists ho > such that, for every < h < ho, 



lim sup — log 

and, more precisely, 



P^(n G r|0 G r) 



EmeN mK h {m) 



-F(0,h) (6.9) 



P^(n G r|0 e r) ^ n ~°° Q A ) H ] * , , 19 (6.10) 

X] meN mK h (m) 4[smh(V2)] 2 1 1 

with Q{.) such that Ej=i QU) n ~°° L{n) / (an a ) . 



It is important to emphasize that, even under assumption (16.6|) . this result would be 
false without the restriction of h small. 

In the light of (16. 9h . it is quite natural to expect (and in some case this can be proven, 
see Section T6.3p that in presence of disorder £ is still proportional to the inverse of the free 
energy, at least close to the critical point. But then, what about ^ av l 

6.2. /j, versus F. To answer this question, we abandon for a while the correlation length 
and we discuss the relation between free energy and another quantity which, due to lack 
of a standard name, we will call simply [i. This was first introduced, to my knowledge, in 
[3] in the context of random heteropolymers: 



fi(P,h) = - lim ^-logE 



1 



ZN,w(P,h) 



(6.11) 



Existence of the limit in our context is easily proven by super-additivity of log Zn )U ((3, h) 
(see |25t Th. 2.5]). An argument similar to (I2.12p gives immediately fj, > while a simple 
application of Jensen's inequality shows that fj,((3, h) < F(f3, h). However, much more than 
this is true: 
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Theorem 6.2. [IT] For every (3 > there exists < c^{f3) , c^{(3) < oo such that 

< <%W ^F$ 2 h) < M> h ) < F (P, h) (6.12) 
if0<h-h c ((3)<c 4 ((3). 

In particular, the bounds in (|6.12|) show that also /i vanishes continuously at the critical 
point, like the free energy. If we call tjf and 77^ the critical exponents associated to the 
vanishing of F and fi for h — * h c {(3) + , Theorem 16.21 implies the following bounds: 

(2 <)vf <Tfr<VF + l, (6.13) 

the inequality in parentheses being valid for (5 > thanks to Theorem 15.81 Just to give a 
flavor of why [i is relevant in the description of the system let us cite the following result. 
Define first of all A at as the largest gap between points of r in the system of length ./V: 

A N := max {\i - j\ : i G r, j € r, {i + 1, . . . , j - 1} n r = 0}. (6.14) 

l<i<j<N 

Then, 

Theorem 6.3. [25J Let (/?, h) £ C. For every e > 0, 

iv lim oo (w^j * ^ - w^yj = 1 m p robamt y- ( 6 - 15 ) 

6.3. Correlation lengths and free energy. To my knowledge, the only case where £ 
and £ av can be fully characterized even in presence of disorder is the one where K{.) is the 
law of the first return to zero of the one-dimensional SRW conditioned to be non-negative. 
In other words, let {Sn}n,=o,i,. .. be the SRW on Z started at 5 = and let P SRW (.) denote 
its law. We define K SRW > + \n) := P 5W (inf{A; > : S k = 0} = 2n\Si > Vi). Go back to 
Section T2.6I for a motivation of this example as a model of wetting of a (1 + l)-dimensional 
substrate. In this case, one has the following 

Theorem 6.4. [41J Let K{.) = K SRW > + {.) and £eZ. For every (3 > and h > h c ((3), 
and, P( du))-a.s., 



7^ = - ^logECJ£ + k,k)= n{(3, h) (6.16) 



l = - lim ^logC^ + fc.A;) =F(0,h). (6.17) 

With respect to Theorem 16.11 this result is much less sharp in that it catches only 
the exponential behavior of the two-point function. However, note that in Theorem 16.41 
h—h c {(3) is not required to be small as in Theorem l6.11 Note also that in Eqs. (|6.16p . (|6.17p 
we have not taken the absolute value of C u (t + k,k): this is because, in this particular 
case, one can prove that this quantity is non-negative [IT]. Finally observe that, in view 
of (|6.12p . the two correlation lengths are different. It would be extremely interesting to 
know whether the two associated critical exponents T]F,r]n coincide or not. 

Remark 6.5. Theorem 16.41 does not coincide exactly with |41| Th. 3.5], e.g., because in 
the latter PS)^ W (.) is the infinite- volume Gibbs measure obtained from the system defined 
in {1, ... , N} letting N — > 00 (cf. footnoted]). However, the proof of [H] extends without 
difficulties to the result we stated above. We remark also that the theorem holds as well 
in the case where K(n) = K SRW (n) := P sw (inf{£; > : S k = 0} = 2n), i.e., the law 
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of the first return to zero of the unconditioned SRW. This follows from the discussion in 
Section EH and from the fact that K SRW (n) = 2K SRW ' + (n). 

Proof of Theorem \6.4\ (sketch). The proof of Theorem 16,41 is based on a coupling 
argument. For simplicity let P + (.) := P SRW (.\S > 0). One can then rewrite the two- 
point function f|6.2[) as 

C u (k,£) = Jim 1 2 

Af^oo Z NtUJ {/3, h) 2 

+i82 rE^ /2 (^n+/ l )(l ( 1} +1 ( 2) )/ _ WD =0 

v {4 =°> i s k =°}J e 

where SW,S® are independent with law P . Since the SRW conditioned to be non- 
negative is a Markov chain, the expectation in the right-hand side clearly vanishes if we 
condition on the event that there exists £ < i < k such that = . But (and here we 
use explicitly the condition Si > and that two SRW trajectories which cross each other 
do necessarily intersect), if the complementary event happens then either or has 
no zeros in the interval {£ + 1, . . . , k — 1}. As a consequence, one obtains 

ICufoO) <2EP^(rn{l,..., h -1} = 0) (6.18) 

and it is not difficult to deduce from (|6.1ip that this probability vanishes like exp(— kfi(/3, h)) 
for k — > oo. For the opposite bound and for the proof of (|6.17p we refer to |41j . 

In the general case where P is not necessarily the law of the returns of the SRW (or, 
in general, of any Markov chain), the available results on correlation lengths in presence 
of disorder are much less sharp and, above all, only correlation length upper bounds are 
known. At present, the best one can prove in general about average correlation length is 
the following: 

Theorem 6.6. [42] Let e > and ((3,h) £ C There exists C x := C x {e,f3,h) > such 
that, for every k E N, 

E\C u (£ + k,£)\ < c (p h)1/Cl exp(-fcCi/x(AM 1+e )- (6-19) 

The constant Ci(e,(3,h) does not vanish at the critical line: for every bounded subset 
B C C one has infV^ C\(e, (3, h) > Ci(B, e) > 0. 

Remark 6.7. The necessity of introducing e > (i.e., of weakening the upper bound 
with respect to the expected one) is probably of technical nature, as appears from the fact 
that for (3 = Theorem 16.61 does not reproduce the sharp results (|6.9p which hold for the 
homogeneous case. 

Observe that Theorem 16.61 is more than just an upper bound on £ av . Indeed, thanks to 
the bound on the prefactor in front of the exponential, Eq. (|6.19p says that the exponential 
decay, with rate at least of order fi 1+e , starts as soon as k 3> t Lle \ 1°S A* I - This observation 
reinforces the meaning of Eq. (|6,19p as an upper bound of order /z^ 1 on the correlation 
length of disorder-averaged correlations functions. 

About the typical correlation length the following can be proven: 

Theorem 6.8. [42] Let e > and (/3, h) € C. One has for every 

MA:, 0)1 < C 2 {Lu)exp(-kC 1 F{(3,h) 1+e ) , (6.20) 
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where C\ is as in Theorem 1 6. 61 while 6*2(0;) := C^w, e, /3, /i) is an almost surely finite 
random variable. 

The proof of Theorems 16.61 and 16.81 relies on a rather involved coupling/comparison 
argument. In simple (and imprecise) words, one first approximates K{.) with a new law 
K{.) which is the law of the returns to zero of a Markov process with continuous trajectories 
(defined in terms of a Bessel process), and at that point the coupling argument of last 
section can be applied. We refer to [32] for full details. 



Appendix A. Two Tauberian results 

For completeness, we include without proof two Tauberian theorems (i.e., results about 
the relation between the asymptotic behavior of a function and of its Laplace transform) 
which we used in Section 15.51 Given a function Q : N — > R, we define for s S R 

Q(s) :=]Te-" s Q(n) 

nSN 

whenever the sum converges. 

We begin with a (quite intuitive) fact: 

Theorem A.l. Proposition 1.5.8] If £(.) is slowly varying and 7 > —1 then 

N ATT+ 1 

Next we state Karamata's Tauberian theorem [8j Th. 1.7.1] which for our purposes may 
be formulated as follows: 

Theorem A. 2. Assume that Q(n) > for every n £ N, that £(.) is slowly varying and 
that p > 0. The following are equivalent: 

Q(s) ~ — ^— ( A - 2 ) 

and 

£q(„)^ W (A.3) 

n=l ^ 

Recall that the function T{z) can be defined, for z > 0, as 



oo 



r(*) = / t*- A e-*cit. 
Jo 

Finally, a theorem relating the Laplace transform of a law on the half-line to its inte- 
grated tail (cf. [HI Corollary 8.1.7]): 

Theorem A.3. Let X be an integer-valued random variables with law P and Q(n) : = 
P(X = n), £(.) a slowly varying function and < a < 1. The following are equivalent: 

l-Q(s) ~°s a l(l/s) (A.4) 



and 
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